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ABSTRACT 

Accreting neutron stars (NSs) often show coherent modulations during type I X-ray bursts, called burst os- 
cillations. We consider whether a nonradial mode can serve as an explanation for burst oscillations from those 
NSs that are not magnetic. We find that a surface wave in the shallow burning layer changes to a crustal inter- 
face wave as the envelope cools, a new and previously uninvestigated phenomenon. The surface modulations 
decrease dramatically as the mode switches, explaining why burst oscillations often disappear before burst 
cooling ceases. When we include rotational modifications, we find mode frequencies and drifts consistent with 
those observed. The large NS spin (ss 270-620 Hz) needed to make this match implies that accreting NSs are 
spinning at frequencies ps 4 Hz above the burst oscillation. Since the long-term stable asymptotic frequency 
is set by the crustal interface wave, the observed late time frequency drifts are a probe of the composition 
and temperature of NS crusts. We compare our model with the observed drifts and persistent luminosities of 
X-ray burst sources, and find that NSs with a higher average accretion rate show smaller drifts, as we predict. 
Furthermore, the drift sizes are consistent with crusts composed of iron-like nuclei, as expected for the ashes 
of the He-rich bursts that are exhibited by these objects. 

Subject headings: stars: neutron — stars: oscillations — X-rays: bursts — X-rays: stars 



1. INTRODUCTION 

Type I X-ray bursts are the result of unstable nuclear burn- 
ing on th e surface of accreting neutron stars (NSs) (s ee re- 
views bv lBildsted[l998t IStrohmaver & Bildstenll2003l) . trig- 
gered b y the extreme temperature s ensitivity of triple-a re- 
actions (|H^MejLj^_^ iJjQrn|ll^Z5l; |WoQsiey_(^_Taai3 LL9_7_6l 
Maras chi & Cavalierell977HJossll977HLamb & Lambll978l) 
They have rise times of seconds with decay times ranging 
from tens to hundreds of seconds, depending on the composi- 
tion of the burning material. When the NS is actively accret- 
ing the bursts repeat every few hours to days, the timescale to 
accumulate an unstable column of fuel. 

Oscillations are often seen in the burst light curves both be- 
fore and after the burst peak ( Muno et alJ2001l and references 
therein). They have frequencies of 270 - 620 Hz (with one 
case of 45 Hz, iVillarreal & Strohmaveil 120041) and typically 
show positive drifts of sa 1-5 Hz. During the burst rise, the 
frequency and amplitude evolution are consistent with a hot 
spot from the burst ignition s preading over the NS surface 
( IStrohmaver. Zhang. & Swankl 19971) . but the positive drift in 
the decaying tail has not ye t been satisfactorily explained. 
iCumming & Bildsten (2000) explored Strohmayer et al.'s 
(1997) hypothesis that this spin-up is simply angular momen- 
tum conservation as the surfa ce layers expand and c o ntract . 
iHevll ( 120001) and lAbramowicz. Kluzniak. & Lasotal d2001l) 
poin ted out the impo r tance o f general relativistic corrections, 
and [Cumming et al. (2002) eventually concluded that this 
mechanism underpredicts the observed drift sizes. These 
works did not resolve the cause of the surface asymmetry at 
late times, long afte r any hot spots should have spread over 
the su rface (Bildsten 1995; Spitko vskv. Levin. & Ushom irskv 
120021) . The asymptotic frequency is characteristic to a given 
object and is very stable over many observations (within 1 part 



in 10 3 . lMuno etalJ2002l) . This fact, along with burst oscilla- 
tions seen from two accreting millisecond pulsars at their non- 
bursting pulsar frequency (SAX J1808.4-3658, Chakrabarty 
etal. 2003; XTE J1814-338, Strohmayer et al. 2003), have 
led many to conclude that burst oscillations exactly indicate 
the NS spin frequency. 

However, it remains a mystery as to what creates the sur- 
face asymmetry in those accreting NSs that do not show 
pulsations in their persistent emission. In Table 1 we sum- 
marize the burst oscillations seen from 12 non-pulsar NSs. 
Since these objects have weaker magnetic fields than the 
accreting pulsars, their burst oscillations may well be due 
to a different mechanism. This hypothesis is supported by 
many differences between the burst oscillations from pul- 
sars and non-persistently pulsating NSs. The non-pulsars 
only show burst oscillations in short (2 -10 s) bursts (ex- 
cluding superbursts, see Table 1), while the pulsars have 
also show n burst o scillations in longer bursts (in XTE 
J1814-338, Strohmayer et al. 2003). The non-pulsars show 
frequency drifts during burst cooling, often late into the 
burst tail, while the pulsars only show drifts during the 
burst rise with no frequency evolution after the burst peak 
dChakrabartv et alJl2003l IStrohmaver et alJl2003l) . The non- 
pulsars have burst oscillations that are highly sinusoidal while 
the pulsars show slight harmonic content (compare the re- 
sults of Muno, Ozel, & Chakrabarty 2002 with Strohmayer 
et al. 2003). Finally, the pulsed amplitude as a function 
of e nergy is different between the two categories of ob- 
jects (iMuno. Ozel. & Chakrabartvl2003HWatts & Strohmaverl 
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An attractive explanation for the burst oscillations in non- 
magnetic N Ss is that they originate in nonradial oscillations 
JHevll 12004) since this is an obvious way to make large 
scale asymmetries in a liquid. Nonradial oscillati ons on 
burstin g NSs were previously studied bv lMcDermott & Taaml 
(1987), but this w as befo re the discovery of burst oscillations 
( Strohmav er et alj|l996l) . They did not incorporate many im- 
portant physical details such as a fast NS spin and the crustal 
interface mode which are crucial to our arguments. The an- 
gular and radial eigenfunctions that are allowed for such a 
mode are se verely restric ted by the many properties of burst 
oscillations. iHevll (12004) identified that the angular structure 
must be given by an m = 1 buoyant r-mode. His arguments 
for this are as follows. The highly s inusoidal nature of the os- 
cillations ( Mun o. Ozel. & Chakrabartvll2002l) . implies an az- 
imuthal quantum number of m = 1 or m = — 1 for the surface 
asymmetry. A mode with frequency wina frame co-rotating 
with the stellar surface is seen by an inertial observer to have 
a frequency 



Density Depth 



w obs = \mCl-cj\, 



(1) 



where SI is the NS spin. Since the frequency of a surface wave, 
ui, decreases as the star cools, the mode must be traveling ret- 
rograde to the spin (m > 0) for a positive frequency drift in the 
observer's frame. The fast NS spin (which we argue is close to 
the burst oscillation frequencies of « 270-620 Hz) alters the 
latitudinal eigenfunctions, so that the angular eigenfunctions 
are no longer given by spherical har mon ics, and ui is simi- 
larly modified (as we summarize in SI4.2I also see Longuet- 
Higgins 1968; Bildsten, Ushomirsky, & Cutler 1996, here- 
after BUC96; Piro & Bildsten 2004). Due to these effects, 
IHevll ( 120041) concluded that r-modes are ideally suited to be 
burst oscillations because they have m > 0, their latitudinal 
eigenfunctions span a wide region around the equator so that 
they may be easier to observe, and the rotational modifications 
to ui result in small drifts. 

The remaining outstanding problem with Heyl's explana- 
tion was in identifying the correct radial structure of the mode. 
The outer NS during burst cooling has three separate layers: a 
hot bursting layer that was heated during the burning, a cooler 
ocean below, and finally the shallow region of the crust (see 
FigureQ. A natural first guess for the cause of the burst os- 
cillation is a shallow surface wave excited in the hot bursting 
layer and riding the buoyancy at the bursting-layer/ocean in- 
terface. One problem with this explanation is that this mode 
overestimates the observed frequency shifts (as the bursting 
layer cools from 10 9 K to a fewx 10 8 K, the mode frequency 
changes by too much, even once rotational modifications are 
included). In addition, this mode cannot reproduce the ex- 
treme stability of the measured asymptotic frequencies (a sur- 
face wave's frequency will vary depending on the temperature 
and resulting composition that is unique to each burst). 

We solve these difficulties by including both the shallow 
surfac e wave and the crustal interface wave ( Piro & Bildsten 
120051 hereafter PB05). The latter is concentrated in the ocean 
and rides the ocean/crust interface. These two modes have fre- 
quencies close enough that they undergo an avoided crossing 
during the burst cooling. The case we focus on here is one in 
which the energy in the surface wave "adiabatically" changes 
into the crustal interface wave as the burst cools. This switch 
solves the two problems listed above because: (1) switching 
to a crustal interface wave ceases the frequency evolution so 
that the predicted frequency drifts match the observed drifts, 
and (2) the final frequency is stable and characteristic to each 
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FIG. 1. — Diagram of the NS surface layers during the decaying tail of 
an X-ray burst. The surface is composed of three regions, a hot layer of 
recently burned fuel with temperature 7J,, a cooler ocean with temperature 
T c , and the outer crust also with T c . The bursting-layer/ocean interface is 
at a depth determined by where unstable burning ignites during the X-ray 
burst, while the ocean/crust interface occurs at the crystallization depth of the 
ocean material. Each interface provides a buoyancy that is a natural place for 
supporting modes. During the beginning of the burst the region just above the 
bursting-layer/ocean interface is convective, but we ignore this complication 
to simplify our analysis. 



object because it only depends on the crust's properties, which 
do not change from burst to burst. 

1.2. Outline of Paper 

In 5j2]we use simple analytic arguments to describe our idea 
for how the surface wave and crustal interface wave interact 
to reproduce the main features of burst oscillations. In ^3] 
we present time dependent numerical models to simulate the 
NS surface during the cooling following an X-ray burst. This 
model is used in |4]to calculate the eigenfrequency spectrum, 
focusing on the surface wave and crustal interface wave, in- 
cluding rotational modifications expected for a quickly rotat- 
ing NS. We also estimate the time dependent oscillation am- 
plitude, showing that it decreases as the surface wave changes 
into a crustal interface wave. In ^5]we compare the observed 
drifts and persistent luminosities with our analytic mode es- 
timates. We consider other modes in 5j6] and examine the 
prospects for their detection. We then conclude with a sum- 
mary of our results in 50 

2. ANALYTIC SUMMARY AND DEPENDENCE ON THE 
NEUTRON STAR CRUST 

We now outline why the shallow surface wave and crustal 
interface wave are consistent with the properties of burst os- 
cillations. This demonstrates which properties of the NS crust 
are constrained by burst oscillation observations. 

A shallow surface wave in the hot bursting layer (see Figure 
^ has a frequency 



2 u ,2 A P 

u s =gH b k , 



(2) 



where g = GM/R 2 is the surface gravitational acceleration, 
H/, = P/(pg) is the pressure scale height at the base of 
the bursting layer, k 2 = X/R 2 is the transverse wavenumber, 
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Ap/p = l-T c Hb/(TbH c ) is the fractional density contrast at the 
burning depth, and pb (p c ) is the mean molecular weight in the 
bursting layer (ocean). Throughout we assume that the ocean 
and crust have the same temperature and composition. In the 
non-rotating limit A = /(/+ 1), where I is the normal spherical 
harmonic quantum number. For a quickly spinning NS, A is a 
function of O and us, which can be accurately calculate d usi ng 
the "traditional approximation" ( as we summarize in 34.21 or 
see BUC96). In agreement with Heyl (2004) we show that 
the only low order, rotationally modified angular mode that 
reproduces all of the properties of the burst oscillations is a 
buoyant r-mode, which has A ~ 1/9 ~ 0.11 in the quickly 
spinning limit (this mode is often iden tified as the / = 2, m = 1 
r-mode from the slowly rotating limit, Piro & Bildsten 2004). 
Using the scalings from equation (0, 

2tt \A b \ 10 9 K / V R 



x 



A 

an 



1/2 



Tc Pb_ 
U p c 



1/2 



(3) 



where A/, and Z/, are the mass number and charge of the ions 
in the bursting layer, respectively, and we use an ideal gas 
equation of state. The final term from the density disconti- 
nuity is k 0.8 near the burst peak. This frequency is consis- 
tent w ith the low-order modes found by McDermott & Taam 
( 119871) . but with A = 2, as appropriate for the / = 1 modes on a 
non-spinning NS that they studied. 

The crustal interface wave has a frequency set by the cool 
NS ocean (PB05), 

cj 2 = p og H c k 2 , (4) 

where po = p/P ~ 10~ 2 (see 34. 1> is the ratio of the shear 
modulus to the pressure at the top of the crust, and the scale 
height, H c , is evaluated at this same depth. Equation is 
the dispersion relation for a shallow surface wave, but with 
a factor of po due to crust compression. For the tempera- 
tures expected at the bottom of the ocean, and using a pressure 
dominated by degenerate, relativistic electrons (PB05) 



2tt \aJ V3x10 8 K 



/10 km 



A 

an 



1/2 



1/2 



where the prefactor is set to match our numerical results, and 
A c is the mass number of nuclei in the crust. 

The shallow surface mode's frequency decreases as the 
layer cools. Once Tb cools to T c « 3 X 10 s K (and using 
fJ'b/^c ~ 1 -0), we find that the surface wave's frequency drifts 
by » 9 Hz, much larger than the observed shifts. This is one 
of the key problems with explaining the burst oscillations with 
only a shallow surface wave. When the crustal interface wave 
is included, it is possible that during the cooling uj s = uj c , at 
which point the surface wave evol ves i nto a crustal interface 
wave (a transition we discuss in 34. 3> . The frequency then 
remains fixed because T c does not change during the X-ray 
burst. The drift then is only Au = uj s - u> c , where u> s is its 
value at the beginning of the burst when 7/, w 10 9 K, result- 
ing in Aui/(2tt) < 5 Hz, much closer to observations. In this 
picture, the size of the drift is approximately set by the com- 
position and temperature of the crust. A crust that is hotter 
or composed of lighter elements has a higher crustal inter- 
face wave frequency, and therefore its drifts are smaller. The 



observed drifts require a crust with T c = (2-6) x 10 8 K for 
A c = 50- 100. We come back to these scalings when we con- 
sider the observed drifts in ^5] 

3. COOLING X-RAY BURST ENVELOPE MODELS 

We numerically calculate the time evolut ion of the cool- 
ing sur face layers and ocean following Cumming & Macbeth! 
(2004). Since the NS radius, R « 10 km, is much greater than 
the pressure scale height H c w (2-4) x 10 4 cm, we approx- 
imate the surface as having a constant gravitational acceler- 
ation, g = GM/R 2 «2x 10 14 cm s- 2 (M/1.4M Q )(10 km/R) 2 
(neglecting general relativity), and plane-parallel geometry. 
We use z and x as our radial and transverse coordinates, re- 
spectively, and in addition find it useful to use the column 
depth, denoted as y (defined by dy = -pdz), giving a pressure 
P = gy. 

The initial profile of the envelope following the burning 
is that of a hot, flux dominated, bursting layer sitting above 
a cooler, lower flux ocean. To mimic this situation, we 
assume a constant flux of Fb « 10 25 erg cm" 2 s" 1 (basically an 
Eddington flux as seen in radius expansion bursts that often 
proceed burst oscillation s) above a bu rning depth column 
of y b = 3 x 10 8 g cm -2 |Bildsten| [l998). The choice of this 
depth mainly affects the cooling timescale, with a deeper yb 
resulting in an extended light curve. This has a small effect 
on the frequencies since for an electron scattering dominated, 

constant flux envelope T/, oc y?' 4 ', so from equation Q, 

1 /8 

u> s oc yJ . Previous studies of constantly accreting NSs have 
shown that the interior thermal balance is set by electron 
captures, neutron e missions, and pycnonuclear reaction in 
the in ner crust (IM iralda-Escude, Paczvnski, & HaenseB 
19901 IZdunik et alj 119921 iBildsten & Brownl Il997t 
Brown & Bildsten! Il998t IBrownl 120001 120041) which re- 



lease « 1 MeV/wip « 10 18 erg g" 1 ( lHaensel & Zdunik1 U990l 
2003). Depending on the accretion rate and thermal structure 
of the crust, this energy will either be conducted into the 
core or released into the ocean such that for an Eddington 
accretion rate up to w 92% of the ener gy is lost to the 
core and exits as neutrinos (IBrownl 12000). " Since we wish 
to investigate how the properties of the crust affect the 
characteristics of the burst oscillations we treat this crustal 
flux, F c as a free variable, and consider its effects around a 
(5) value of 10 21 erg cm" 2 s" 1 (as expected forM » 10~ 9 M, 



© yr . 

about one-tenth the Ed dington rate) . 

We use the results of Schatz et al. ( 200 lj) and lWooslev etafl 
(2004) to set the compositions. Table 2 lists the key properties 
of our models, which sample a range of conditions expected 
for the surfaces of bursting NSs. We consider only one nu- 
clear species per layer, where that species is chosen to repre- 
sent the wide range of elements that actually exist. In He-rich 
bursts, the burning predominantly produces a-elements and 
cooling lasts ~ 10 s as observed with burst oscillations. We 
expect more mixing dur ing these bursts due to convection at 
the ig nition depth ( Woosl ev et all2004l lWeinberg & Bildsten 
120051) . so that fib/ ' Pc ~ F Models 1 and 2 are meant to rep- 
resent this type of burst, with Model 2 exploring the changes 
that occur for a large T c , On the other hand, bursts with a com- 
position closer to solar (mixed H/He bursts) produce a heavy 
ocean of ashes from rp-process burning and also less mixing 
so that pb /fe 5, !■ These bursts occur on objects such as GS 
1826-24 (Gallowav et al. 2004) and have never shown burst 
oscillations. Model 3 represents this case, but without the ex- 
tra heat source of the rp-process burning. This only changes 
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the timescale of the cooling and frequency evolution, and does 
not affect the calculated mode frequencies nor drifts. 

The density of the liquid/solid transition, p c , when crystal- 
lization begins is set by the dimensionless parameter 



r= 



(Zee? 



a c k B T c 



= 169 



3 

7^ 



1/3 



10 9 g crrr 



1/3 



, (6) 



where is Boltzmann's constant, T c $ = r e /10 8 K, Z e is the 
charge of the crust nuclei, and a c = (3/47r«;. c ) 1//3 is the average 
ion spacing with n l c the ion number density in the crust. This 
transition occurs at V w 173 (Farouki & Hamaguchi 1993 and 
references therein), implying a density at the top of the crust, 



p c = 1.1 x 10 g cm 



3 



-V- 

173 / 



(7) 



Since material in the crust is simply ocean material that was 
advected there by accretion, we assume that the ocean and 
crust have the same composition. 

We next construct our envelope models using the parame- 
ters summarized above and follow their evolution forward in 
time. Cooling is described by the heat diffusion equation with 
no source terms, 



dT 

Cp ~dF 



OF 



(8) 



where c p is the heat capacity at constant pressure, and the 



radiative flux is 



F = 



4acT 3 dT 
3k dy ' 



(9) 



where a is the radiation constant and k is the opacity . The 
opacity is set using electron-scattering (Paczvhski 1983), free- 
free (Clayton 1993 with the Gaunt factor of Schatz et al. 
1999), and conductive opacities (Schatz et al. 1999 using 
the basic form of Yakovlev & Urpin 1980), and we assume 
that the crust opacity is set in the same way. We solve fo r 
p using the analytic equation of state from PaczvriskH ( 119831) . 
The initial temperature profile is found by integrating equa- 
tion (|9jl with the flux profile described above. The envelope 
is then evolved in time according to equation l|8} using fi- 
nite differencing techniques. We use a grid that is uniform in 
sinh -1 [log(y/yt,)] (|Cummin2 & Macbeth 2004J, so as to prop- 
erly resolve the temperature jump at the burning depth. Figure 
|2]shows the initial profile along with the profiles at time steps 
from 0.1 - 10.0 s for Model 1. The majority of the evolution 
is concentrated in the outer hot bursting layer with a negligi- 
ble thermal wave diffusing into the ocean and crust, due to the 
long therm al time at these depths (consistent with the work of 
Woos lev"et al.12 004). These general features are shared by all 
of the models we consider. 

4. NONRADIAL MODES DURING X-RAY BURST COOLING 

We now calculate the time dependent spectrum of mode 
frequencies, focusing on the surface and the crustal interface 
waves. Even though the NS is spinning quickly, VL/{2n) « 
270-620 Hz, the buoyancy of the envelope dominates over 




10 7 10 8 10 9 

Column Density (g cm -2 ) 

FIG. 2. — Temperature, density, and flux for a He-rich bursting NS en- 
velope (Model 1). The initial profile is a flux of 10 25 erg s _1 cm -2 above a 
column of 3 X 10 s g cm~ 2 , with a flux of 10 21 erg s cm~ 2 below (thick solid 
lines). The curves showtime steps ofO.l s (solid), 0.3 s (dotted), 1.0 s (short- 
dashed), 3.0 s (long-dashed), and 10.0 s (dot-dashed). The majority of the 
time dependent change occurs in the outer hot bursting layer where the local 
thermal time is short. These profiles extend down into the crust (which begins 
at a column of 5.3 X 10 12 g cm -2 ), but here we just focus on the region where 
the main changes occur. 



Coriolis effects and determines the radial structure of the 
modes (BUC96). This is because N 2 ~ 10 1() s" 2 > RQuj/H c - 
10 6 s~ 2 , where the internal buoyancy of the envelope is mea- 
sured by the Brunt-Vaisala frequency 



N 2 = Jd\ogp 



1 d\ogP 

\ dz T\ dz 



(10) 



and Ti = (<91ogP/<91og p) s is the adiabatic exponent. Further- 
more, since the modes are squeezed into a thin outer layer 
of thickness « H c <C R, they have predominantly transverse 
amplitudes. Together these properties justify accounting for 
the effects of rotation by using the "traditional approxima- 
tion" (BUC96). This simplification allows the angular part 
of the perturbation equations to be separated from the radial 
part, and solved as an eigenvalue equation for A (the "effec- 
tive wavenumber"), where A is a function of Q, and u>. For 
this reason, we focus on the radial equations first, and then 
introduce the angular and spin dependent parts afterwards in 



4. 1 . Radial Mode Equations 

Solving for the radial eigenfunctions requires considering 
the governing mode equations both in the bursting layer and 
ocean, and also in the NS crust. The bursting layer and ocean 
have no shear modulus, so adiabatic perturbations in this re- 
gion are described by the nonradial oscillations equations for 
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an inviscid fluid in hydrostatic balance and plane-parallel ge- 
ometry ( Bi ldsten & Cutlerl 19951 hereafter BC95) 



d& _ & _ (gHk 2 _ 1 \ SP 



dz TiH \ 



d 5P 
~dz~P 



r, 



1 SP 
H~P 



TiJ P ' 
N 2 \ & 



(11) 



(12) 



where is radial displacement and SP/P is the Eulerian pres- 
sure perturbations. The Eulerian perturbations are related to 
the Lagrangian perturbations by AP = SP-!; z pg. When cal- 
culating g-modes (including the shallow surface wave), equa- 
tions (II \\ and (11 21 are sufficient since these modes are ex- 
cluded from the NS crust (BC95). The only exception is the 
crustal interface mode (PB05), where it is crucial to solve the 
mode equations with a nonzero shear modulus, /i, in the NS 
crust (BC95 ). For a classical one-co mponent plasma the shear 
modulus is (IStrohmaver et all!991l) 



0.1194 



n Uc (Z c ef 



(13) 

1 + 0.595(173/T) 2 a c ' 

where T is given by equation (|6}. Since the pressure in the 
crust is dominated by degenerate, relativistic electrons we 
rewrite /x as 

1.4 xlO" 2 /Z,,^ 2/3 



P l+0.595(173/r) 2 V30 



(14) 



In the crust fi/P is fairly independent of temperature (except 
for a small dependence in the factor of T in the denominator), 
so we substitute fio = /i/P and assume that /j,o is constant with 
depth (PB05). 

4.2. Rotational Effects 

Before the mode equations can be integrated, we must con- 
sistently set A with respect to ui and 57. The dispersion relation 
that relates these three variables is described by BUC96 (from 
the work of Longuet-Higgins 1968), and results in a large va- 
riety ofangular eigenfunctions for a given radial structure. 
iHevll ( 120041) showed that very few of these (just the m = 1 
buoyant r-modes) match the properties of burst oscillations. 
We summarize these results here, both for completeness and 
b ecause we consider the presence o f additional modes in 5j6] 

Maniopoulou & Andersson (2004) have extended the tradi- 
tional approximation to include general relativity, which in- 
troduces corrections to the co-rotating mode frequencies, u>, 
due to frame-dragging and red-shifting. Each of these ef- 
fects contribute differently depending on the angular eigen- 
function under consideration, and combined they can decrease 
u> by as much as 20%. Since the buoyant r-modes have fre- 
quencies independent of spin (as we describe below), we ex- 
pect frame-dragging to be negligible and red-shifting to be 
the main change to the normalization of u> (as is fou nd for 
the Kelvin mode by Maniopou lou & Anderssonl2004l) . From 
equations (0 and (|3J we see that such corrections are degen- 
erate with slight changes in the radius and temperature of the 
NS, and therefore we do not include them. 

There are two sets of solutions for A, which are both typi- 
cally presented in the literature as a function of the "spin pa- 
rameter," q = 20,/uj (so that for a quickly spinning star q 3> 1) 
The first set is comprised of the rotationally modified g-modes 
(those with m ^-l) and the Kelvin modes (to = —[), which are 
shown in Figure [5] In the slowly rotating limit, q <C 1 , we 
find A«/(/+l)~2or6 (for I = 1 or 2, respectively), but as 



100 - 



10 



0.1 




q = 2fi / cj 



FIG. 3. — Effective wavenumber A as a function of the spin parameter 
q = 2Q/w for the rotationally modified g-modes (m^—1) and Kelvin modes 
(m = -/). Solid lines denote modes for which / = 1 in the non-rotating limit 
with m = — 1, 0, and 1 (from bottom to top). The dashed lines denote modes 
for which / = 2 in the nonrotating limit with m = —2, —1, 0, 1, and 2 (from 
bottom to top). For q 3> 1, A oc q 1 , except if m = —I, which asymptotes to 
A = m 2 (the Kelvin modes). Unfortunately, none of these dispersion relations 
are consistent with the burst oscillations because they either have frequency 
shifts much too large (in the case of the rotationally modified g-modes) or 
shifts with the wrong sign (in the case of the Kelvin modes). 



the spin increases A splits depending on its to value. It then 
asymptotes to scaling as A oc q 2 (in the case of the rotation- 
ally modified g-modes) or as A = m 2 (in the case of the Kelvin 
modes). 

We present these solutions to emphasize that none of these 
modes are consistent with burst oscillations. The rotationally 
modified g-modes all have A > 1 when the star is rotating 
quickly. This implies large mode frequencies (see eqs. and 
(5|). and frequency drifts larger than observed. (In <j6] we ex- 
plore whether these modes exist in burst cooling light curves, 
but are not seen because of their large frequency drifts.) The 
Kelvin modes are inconsistent because they are all prograde 
(to < 0), so that as the surface cools an observer would see a 
decreasing frequency. 

The other set of solutions are a group of modes unique to the 
case of a rotating star, the r-modes, which we plot in Figure 
13 These occur as zero frequency solutions for a non-rotating 
star, corresponding to incompressible toroidal displacements 
on the stellar surface. When the star is rotating, Coriolis ef- 
fects turn these solutions into normal modes of oscillation. In 
the slowly rotating limit the NS spin and mode frequency are 
directly proportional, 

7.m£l 

(!Saioll982tfLee & Saioll986T) . which are the well-studied in- 
ertial r-modes. More relevant for our work here is the case 
when q 1, which are the buoyant r-modes. For q 3> 1 the 
r-modes with to = / exhibit A oc q 2 , similar to the rotationally 
modified g-modes from Figure|3] and thus also give frequency 
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FIG. 4. — Effective wavenumber A as a function of the spin parameter 
q = 2U/lu for the r-modes. Solid lines denote modes for which m= 1, with 
values of / (identified in the slowly rotating limit using eq. 1151 ) of 1, 2, and 
3 {from left to right). Dashed lines denote modes for which m = 2, with 1 = 1, 
2, and 3 (from left to right). We favor the / = 2, m = 1 mode (thick line) as 
explaining the burst oscillations, which has A m 1/9 ~ 0.11 in the quickly 
rotating limit. 



shifts too large for burst oscillations. More promising are the 
r-modes with m < I, which have A < 1 and therefore small 
frequency shifts. Furthermore, since all r-modes have m > 0, 
they are traveling retrograde with respect to the NS spin, and 
provide increasing shifts as the NS cools. We favor the 1 = 2, 
m = 1 mode (denoted with a thick line in Figure |4} over any 
of the other A< 1 r-modes because: (1) m = 1 is implied 
from the observations and this is the lowest order m = 1 mode 
that gives frequencies and shifts consistent with burst oscil- 
lations and (2) any higher order A < 1 r-mode has multiple 
bands of hot and cold regions at different latitudes, so that 
it should be more difficult to observe than the I = 2, m = 1 
mode. It is rather remarkable that with the number of rota- 
tionally modified modes that exist, this is the only one that 
fits all the required properties ! For this mode A « 1 /9 « 0. 1 1 
in the quickly rotatinglimit, as we used when estimating the 
mode frequencies in |5] 

A further reason for favoring the buoyant r-modes i s that 
their latitudinal eigenfunctions may be easier to observe 
2004). Modes with A 1 have their eigenfunction squeezed 
near the equator within an angle cos 6 < l/q, where 9 is mea- 
sured from the pole. On the other hand, the m < I r-modes 
and Kelvin modes have much wider eigenfunctions that span 
most of the surface. This means that given a fixed equato- 
rial perturbation, the pulsed fraction is larger for these latter 
modes. 

4.3. Eigenfunctions and Frequencies in Rotating Frame 

Normal modes of oscillation are found by assuming AP = 
at the top boundary, which is set at a depth where the local 
thermal time is equal to the mode period (f t h = 2tt/u), where 
f t h = c p yT/F). This top condition, though not unique, is 
fairly robust since little mode energy resides in the low den- 
sity upper altitudes (BC95). This is especially true for the 



two modes we study here, which have their energy concen- 
trated at their respective interface. We numerically integrate 
the mode equations, shooting for the condition that £ z « t; x rs 
deep within the crust, where £ A is the transverse displace- 
ment. This may not be the case in a more realistic calcu- 
lation of the crust, but as long as we set this bottom depth 
deep enough, w 10 17 g cm" 2 , we recover the asymptotic so- 
lutions of PB05. When we change from integrating the non- 
viscous mode equations to those appropriate in the crust, we 
must properly set boundary conditions (PB05). 

When solving the mode equation s we use an effective wave 
number A « 0.11, as described in 34.21 There are many so- 
lutions to the mode equations, which form a complete basis 
set of functions. We focus on the low frequency solutions (as 
opposed to the /-mode or /?-modes) because they are a better 
match to the observed frequencies. These solutions are or- 
dered from high to low frequency, with each successive mode 
having an additional node in its radial eigenfunctions. In Fig- 
ure [5] we plot the radial and transverse displacement eigen- 
functions (£, and £ v , respectively, given in cm) for the first 
three low frequency solutions from Model 1 after 0.1 s of 
cooling. In addition we plot the energy per logarithm column 
density, 



dE 



UirRWey, 



« log J 

where £ 2 = |£t| 2 + |£z| 2 is the total displacement. This indi- 
cates where the kinetic energy of the mode is concentrated. 
We normalize the total integrated energy of each mode to be 
5 x 10 36 ergs, 1 0~ 3 of the total energy released by unstable 
nuclear burning (Bildstenl ll998l) . At 9.14 Hz (where this is 
the frequency in the rotating frame on the NS surface, lo) we 
find a mode with a single node in its eigenfunctions (shown as 
a cusp in £ z and less apparent in £ A because there is always a 
node at the crust at a column of y c = 5.3 x 10 12 g cm -2 ). Since 
the energy is concentrated in the bursting layer this is identi- 
fied as the shallow surface wave. The next mode at 4.30 Hz 
has an additional node and therefore a different parity than the 
surface wave. Its energy is distributed very differently, largely 
concentrated at the bottom of the ocean, so it is identified as 
the crustal interface mode. All modes with lower frequencies 
are g-modes, and we show the first g-mode (at 3.10 Hz) with 
three nodes. 1 These modes are trapped in either the burst- 
ing layer or ocean and are identified by a relatively constant 
energy distribution in each region (the example shown is the 
prior case). 

Since both the surface wave and the first g-mode are con- 
centrated in the bursting layer, they each exhibit frequency 
drifts as the bursting layer cools, making them both attractive 
as being burst oscillations. We therefore consider the radia- 
tive damping time for each mode to narrow down the choice 
of radial eigenfunction. The rate of energy loss averaged over 
one oscillation cycle from an adiabatic p erturbation is esti - 
mated from the "work integral" jCoxll980l:IUnno et alJ 19891) . 
which, when there is no energy source, is 

dE uj f , /T 3 -l AP A / dF\ „ , 

where Tj,— 1 = {d\ogT / d\ogp) s . The integral is negative 
when a mode is damped and positive when a mode is excited. 

1 These g- mode s should not be confused with the rotationally modified g- 
modes from 34 21 Here we are only considering the radial eigenfunction. In 
the previous case we were discussing the angular eigenfunctions, which can 
be applied to any radial eigenfunction. 
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FIG. 5. — Displacement eigenfunctions (given in cm) and energy density (measured in ergs and given by eq. Ildl ) of the first three low frequency solutions. 
These are from Model 1 after 0.1 seconds of cooling and all with A = 0.11. The frequencies are in the rotating frame on the NS surface, u>, From high to low 
frequency (left to right), these are the shallow surface wave, the crustal interface wave, and the first g-mode. All lower frequency solutions are also g-modes. 
Each lower frequency mode contains an additional node, which are shown by the cusps because we plot the absolute value of the displacements. The transverse 
displacement, £. v , always has a node at the crust at a column of 5.3 X 10 12 g cirT 2 . 



From this and the total mode energy (integrating eq. 1161 ) we 
can estimate the damping e-folding time of a mode's ampli- 
tude, 



fdamp — 



dE 



-dlogy 





dE 


-l 


) 


dt 





(18) 



dlogy 

where the absolute value has been included to make this 
timescale positive (since all modes we consider are damped). 

Figure |6] shows fd am p for both the surface wave and the first 
g-mode as functions of time since the burst peak. To ac- 
curately calculate the integral from equation d!7l > we must 
drop derivatives of AP/P near the surface because our im- 
posed boundary condition of AP = results in unphysical lo- 
cal mode excitation (for a d escription of ano ther way to set 
this boundary condition see iGoldreich & Wull 19991) . The g- 
mode is damped on a timescale an order of magnitude faster 
than the surface wave. Though both the first g-mode and the 
surface wave have energies concentrated in the bursting layer 
(see far left and right panels at the bottom of Figure each 
mode's energy distribution is different within this region. Ra- 
diative damping predominantly takes place at the top of the 
atmosphere (Piro & Bildsten 2004), so that a mode with more 
of its energy distributed there is more damped. It is some- 
what difficult to see because of the large dynamic range we 
show, but the energy of the g-mode is approximately constant 



in the bursting layer, while the shallow surface wave's energy 
varies by almost an order of magnitude in the bursting layer 
(compare the bottom panels on the far left and right of Fig- 
ure |5}. The ratio of damping timescales is therefore approx- 
imately the ratio of energies between the top and bottom of 
the bursting layer. After w 0.5 s the g-mode has a damping 
time shorter than the time for the burst to cool, so that the am- 
plitude of this mode is exponentially damped. Higher-order 
g-modes have different damping times depending on where 
they are trapped. Those trapped in the bursting layer (like the 
example shown) all have similar, short damping times. The g- 
modes that are trapped in the ocean have much longer damp- 
ing times (fdamp ~ 5— 10 s), but suffer from very small surface 
amplitudes. The crustal wave has an even longer damping 
time (fdamp > 10 3 s), but also has the problem of a small sur- 
face amplitude. We therefore favor the surface wave as the 
radial eigenfunction for the burst oscillations because of its 
combined attributes of having a long damping time coupled 
with a large surface amplitude from a small input of energy. 

We next consider how the mode frequencies evolve with 
time since burst peak. In Figure0we plot the first two eigen- 
frequencies as a function of time for Model 1. The mode 
parity does not change along each continuous line, with one 
node present for the upper line and two nodes for the lower 
line. Initially the upper frequency is the surface wave and 
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FIG. 6. — Damping time, f(i amp , as a function of the time since the burst 
peak, for the mode with a single node (which is initially a surface wave) and 
the mode with three nodes (the first g-mode), both from Model 1. The dashed 
line shows were the damping time is equal to the time since burst peak. After 
only rj 0.5 s the g-mode is exponentially damped on a timescale shorter than 
0.5 s, so that it will not be observable. 



the lower is the crustal interface wave, consistent with equa- 
tions and 0. The surface wave's frequency decreases as 
the bursting layer cools, eventually running into the crustal 
mode. This creates an avoided crossing between the two 
modes (which we highlight in the inset panel of Figure 0. 
The eigenfunctions therefore change considerably along each 
line from being characteristic of a surface wave to a crustal 
interface wave (and vice versa). As long as the evolution of 
the envelope is sufficiently slow, once excited a mode will 
most likely evolve along one of these branches and not be 
able to skip to the other because of differences in parity. We 
do not show this rigorously, but think this is a reasonable 
conclusion by considering analogies with perturbation the- 
ory problems in non-relativistic quantum mechanics. In cases 
where the background changes quickly, on a timescale shorter 
than ~2tt/w, mixing is highly possible as determined by tak- 
ing inner products between initial and final states on either 
side of the avoided crossing. On the other hand, in the case 
we consider here the background changes very slowly, as is 
supported by the fact that d(l/u))/dt < 10" 2 <C 1 (which we 
calculated to check our conclusions). This implies that the 
change can be treated as adiabatic, and the modes are not ex- 
pected to mix at the crossing. 

In Figure[8]we plot the energy of the single-node wave for 
a sample of time steps from Model 1 . We normalize the total 
integrated energy of each mode to be 5 x 10 36 ergs. Initially, at 
4.0 s, the mode's energy is concentrated in the recently burned 
portion of the upper envelope. As time progresses the energy 
goes deeper into the star, so that by 12.0 s, when the burst is 
nearly over, the energy is concentrated at the bottom of the 
ocean, just as expected for a crustal mode. 

4.4. Observed Frequencies 

We now focus on the evolving mode with a single node 
because it is initially a surface wave, which as discussed in 
34.31 is the mode most least likely to be damped and also 
have an appreciable amplitude. In Figure |9] we use equa- 
tion Q to plot the frequencies seen by an inertial observer 




Time (s) 

FIG. 7. — The frequency evolution of the surface wave and crustal interface 
wave using Model 1. The upper curve follows the mode with one node in 
its eigenfunctions, while the lower curve follows the mode with two. The 
qualitative shapes of the eigenfunctions vary considerably along each curve. 
Initially, the higher frequency mode corresponds to a shallow surface wave, 
while the lower frequency mode is a crustal interface wave. As the lines move 
to the right toward later times, the surface wave is sensitive to the cooling sur- 
face layers and decreases, while the crustal wave does not evolve. Eventually 
the two modes meet at an avoided crossing, which we display in an inset 
panel. Past the avoided crossing, although the parity has not changed, the 
higher frequency mode is now like a crustal wave, while the low frequency 
mode is a surface wave. 



(solid lines) for a NS spinning at 400 Hz (dotted line) for each 
of our three models. We also show the light curves with a 
dashed line. In the top panel we plot the results from Model 

1 , This shows that the NS spin is rs 4 Hz above the burst os- 
cillation, just the frequency of the crustal mode. In the next 
two panels we plot the frequency evolution of our other two 
models, along with their respective light curves. In Model 

2, the flux coming through the crust is a factor of 3 higher, 
providing a significantly hotter crust. This causes the crystal- 
lization to happen at much larger depths, resulting in a larger 
crustal wave frequency. The mode crossing therefore happens 
earlier, which can be seen by the smaller frequency drift, as 
expected from the analytic considerations in ^2] In Model 

3, there is a slightly larger discontinuity at the burning depth 
due to the heavy ashes from rp-process burning as expected in 
H/He mixed bursts. This makes the evolution of the surface 
wave much slower, so that in this case it does not transition 
into a crustal interface wave, even after 14 seconds of cool- 
ing. This shows that an avoided mode crossing is not a robust 
property exhibited by every model. Observations of burst os- 
cillations from such a NS would show more variance in its 
asymptotic frequency because of this feature. On the other 
hand, since Model 3 represents an rp-process burst and such 
bursts do not exhibit burst oscillations, this has not been ob- 
served. It is not clear from our models why the He models that 
show mode crossings are the ones observed, while our H/He 
models that don't show crossings are not, especially since fre- 
quency drifts like that shown in the bottom panel would be 
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FIG. 8. — Energy of the single-node mode per unit logarithmic pressure 
at different time steps, from Model 1 . Each mode is normalized to a total 
energy of 5 X 10 36 ergs, 1(T 3 of the total energy released in an X-ray burst. 
At early times the mode looks like a shallow surface wave, with its energy 
concentrated near in the bursting layer, but as it passes through the avoided 
crossing (see Figure 171 its energy moves continually deeper and becomes 
more like a crustal interface wave. The sharp feature at a column of 3 X 
10 8 g cnT 2 is due to the density discontinuity at the burning depth. Below a 
column of 5.3 X 10 12 g cnT 2 the crust begins and the mode's energy drops 
off quickly. 
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observable if they existed. We discuss these issues further 
when we conclude in $7\ 

As noted by Heyl (2004), since A « 0.1 1 as long as q > 10, 
any spin NS will exhibit similar frequency evolution as long 
as f2/(27r) > 5u) s /(2ir) w 50 Hz. This makes it easy to com- 
pare models with different spins, and constrains the explana- 
tion of burst oscillations as modes since it requires that the 
frequency shifts be spin independent. T his is counter to the 
studies of Cummins & Bildsten (2000), which predict that 
the spin and drift size should be correlated. Although EXO 
0748-676 is near this critical frequency (see Table 1), its drift 
(assuming that it is not yet a crustal interface wave) is still not 
heavily affected. We find a ss 4 Hz drift instead of a 5 Hz drift 
at higher rotation rates, so we still predict a spin frequency of 
* 49 Hz for this object. 

4.5. Surface Amplitude Evolution 

As the surface wave evolves into a crustal wave, we ex- 
pect the surface amplitude to decrease dramatically since the 
crustal wave's energy is mostly at the bottom of the ocean. 
We estimate SP/P at our "surface" where the modes are no 
longer adiabatic, t± = 2tt/uj, denoted (SP/P) S , which provides 
a general idea for how the amplitude is changing. BC95 show 
that the Eulerian perturbation SP/P is constant in the non- 
adiabatic region, so that we only need to calculate it at our top 



FIG. 9. — The observed frequency evolution of the single-node mode on 
Models 1, 2, and 3 (solid curves). We use a NS spin of 400 Hz (dotted line). 
For Model 1, this shows that while the mode is a surface wave it evolves 
upward in frequency, but once it becomes a crustal wave it is very stable at a 
frequency of u) c /(2n) rj4Hz below the NS spin. In Model 2 the frequency 
shift ends very early and is only 3 Hz because the crust is much hotter due 
to a larger crustal flux. On the other hand, in Model 3 we find that an avoided 
crossing between the surface wave and crustal wave does not occur while the 
luminosity is at an observable level, so that the frequency never becomes flat. 
Observations of a NS with oscillations like this would show more variance 
in their asymptotic frequency in comparison to the other models we consider. 
The dashed line is the surface luminosity as a function of time, assuming a 
radius of 10 km. 



boundary to estimate its value at the photosphere. Approxi- 
mately 5 x 10 39 ergs of energy is released in an X-ray burst, 
so we assume some small fraction of this is able to power 
the mode. Since we do not know how much energy actually 
leaks into the mode, we simply study values which result in 
reasonably sized perturbations. A key conclusion is that only 
a small fraction of the available energy is needed, ~ 10~ 3 , to 
get perturbations of order unity. We assume that this energy is 
conserved throughout the evolution of the mode, an approxi- 
mation that is justified because the damping time of the mode 
is always longer than the time that the burst has been cooling. 
In Figure^] we plot the surface perturbation as a function 
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FIG. 10. — Evolution of the surface pressure perturbation for the single- 
node mode from each of our models. In each case we assume a certain frac- 
tion of the X-ray burst's 5 X 10 3 ' ergs of energy has gone into the modes, and 
is then conserved in the mode as the envelope cools. The fractions we use 
are 10~ 3 (model 1; solid line), 4 X 10~ 3 (model 2; dashed line), and 2 X 1CT 4 
(model 3; dotted line). This demonstrates that appreciable amplitudes can 
be expected, even if only a small fraction of the burst's energy goes into the 
modes. This also highlights how the amplitude quickly dies off once the 
shallow surface wave turns into a crustal interface mode. 



of time for all three models. In both Models 1 and 2 the per- 
turbation drops off dramatically once the surface wave turns 
into a crustal interface wave, especially in the case of Model 
2, which has a thicker ocean. If the top of the crust is too 
deep the resulting surface amplitude will be too small to ex- 
plain burst oscillations. The total displacement is about con- 
stant down to the crust with £ ss £ x w £ Z R / h since the modes 
are nearly incompressible. From our top boundary condi- 
tion of AP = we find (&/A)i = (5P/P) S , so that together 
£ ~ (SP/P) S R, Using the scalings from equation d!6i . along 
with the crustal mode frequency from equation (|5} and the 
density of crystallization from equation 0, the total crustal 
wave energy at late times is approximately 



£totai, c ~ 5 x 10 36 ergs 
A 

x " 



0.11 



64 

A- 
(SP/P) 



30\ 
Z~c) 



20/3 



R 



10 km 



0.1 



3 x 10* K 



, (19) 



where the prefactor has been set to match the numerical re- 
sults. The predominant effect is that the surface amplitude 
is very sensitive to the depth of the crust so that (SP/P) S oc 

20 /3 — 5/2 

Z c T c for fixed energy. In Model 3 the switch to a 
crustal interface wave never occurs, so that the amplitude 
always stays large. An interesting feature of all the mod- 
els is that at early times the amplitude actually increases. 
This is because as the frequency initially decreases the dis- 
placements must increase to maintain the same energy. Sug- 
gestively, the time evolutions of the amplitudes, especially 
for Model 1, are similar to the amplitudes measured by 
iMuno. Ozel. & Chakrabartvl Eool . 



4.6. Magnetic Field Limits 

The lack of persistent pulsations from the non-pulsar burst 
oscillation sources implies a smaller magnetic field for these 
objects. To prevent channel ed accretion the surfa ce dipole 
field must be less than jFrank. King. & Rainell99l 



B acc = 5.1 x 10 7 G 
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,lO-M yr-<< ' (20) 

Both the shallow surface wave and the crustal interface wave 
exhibit large shears in their transverse displacements. Fol- 
lowing BC95, the maximum magnetic field before the waves 
would be dynamically affected is B 2 « 87t// 2 /xj 2 , where H is 
chosen as appropriate for each mode. For the surface wave 
this gives a limit (using eq. |3*|). 



B s = 8.8 x 10' G 



yb 



2Zt, Tb 

1b~l¥~Kj \3x 10 8 gcm" 
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1.4M Q A 
M 0.11 
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1/2 



,(21) 



surprisingly close to B acc considering that these two limits 
result from vastly different physical mechanisms. For the 
crustal interface mode (using eqs. and Q), 



B c = 1.2 x 10 1U G 



64 
A~c 
1.4M 



11/3 
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(22) 



where we assume T = 173. This is much higher than ei- 
ther of the two previous limits because the crustal inter- 
face wave "lives" at such high pressures. The lack of 
persistent pulses from the majority of X-ray bursters is 
likely because of a weak field in the burs ting layer and 
ocean ( Cumming, Zweibel, & Bildsten 2001), making our 
non-magnetic mode calculation adequate for now. 

5. COMPARISONS WITH OBSERVED DRIFTS 

Brown, Bildst en. & Rufledgd dl998l) showed that the core 
of a transiently accreting NS reaches a steady-state set by nu- 
clear reactions deep in the crust, so that the luminosity de- 
puting the crust is given by the time-averaged accretion rate 
(Brown 2000). This means that crustal temperatures can be 
inferred from an observable property, namely the average per- 
sistent luminosity. Given the discussion from |2j we expect 
objects with small observed drifts to have larger persistent lu- 
minosities. 

In Figure ^] we plot predicted frequency drifts using 
the difference of equations and 0. We relate T c to 
the persistent luminosity, L pels , using the crustal models of 
Brown (2005, private communication) for compositions of 
56 Fe (thick solid line) and 106 Pd (thick dashed line). These 
m odels ar e base d on the calculational techniques described 
in iBrownl J2004I) . We then plot the observed luminosities 
and largest observed frequency drifts from Table 1 . The to- 
tal range of observed luminosities is shown by the horizon- 
tal bars with a point placed at the average luminosity. Open 
circles indicate objects for which burst oscillations are coin- 
cident with radius expansion bursts, while solid triangles in- 
dicate those that are not. We use an Eddington luminosity of 
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FIG. 11. — The frequency drifts, Aaj/(27r), as a function of the persis- 
tent luminosity, Lpets (in units of L^a = 2.5 X 10 38 erg s -1 ). Two different 
crustal models are plotted (Brown 2005, private communication) with com- 
positions of 56 Fe (thick solid line) and 10fl Pd (thick dashed line). We assume 
that Aj/Zj = 2, [ib/fic - 1-0, and that initially Tj, = 10' K. We compare this 
with the largest observed frequency drifts of the non-pulsar, bursting NSs (see 
Table 1). The observed drifts should be viewed as lower limits since observa- 
tional or excitation-related effects may result in the initially observed surface 
wave having T\, < 10 9 K. The horizontal bars indicate the range of luminosi- 
ties shown by each object, with a point at the average luminosity. Open circles 
indicate NSs that show radius expansion bursts at the same time as burst os- 
cillations, while solid triangles indicate NSs that do not. 4U 1728 — 34, Aql 
X— 1, and 4U 1636 — 536 have shown both cases. 



Lf m = 2.5 x 10 38 erg s -1 so as to be consistent with the work 
of iFord et alJ (12000). The observed drifts should be viewed 
as lower limits since the surface wave may not be excited un- 
til the bursting layer's temperature is below the 7j = 10 9 K 
that we use here. Qualitatively, the higher luminosity NSs 
show smaller drifts, as expected from our models. Further- 
more, lighter crustal material is favored for explaining the 
observations. This is consistent with the cooling timescale 
of the bursts (see Table 1), which suggest He-rich bursts that 
produce iron-like nuclei ( and would not produce rp-process 
elements such as 106 Pd). In general, these objects are not all 
expected to have the same crustal composition, and this can 
be studied by making further comparisons with detailed mod- 
eling and consideration of each system's bursting and binary 
properties. 

6. COULD OTHER ANGULAR MODES BE PRESENT? 

In 34.21 we considered the multitude of rotationally modi- 
fied solutions available, all with the same single node radial 
structure, and identified only one that satisfied the properties 
of burst oscillations. The other low angular-order modes re- 
sulted in frequency shifts that were too large, and sometimes 
in the wrong direction. Nevertheless, these modes may be 
compelling to study since a priori we do not have a reason 
to exclude them, and more sensitive observations may re- 
veal these additional modes. Such a measurement of multi- 
ple modes would revolutionize our knowledge of an accreting 



NS's mass, radius, spin, and surface layers. 

We calculate the time dependent frequencies from Model 
1, this time considering a variety of the low angular-order 
modes, all with the same single-node radial eigenfunction as 
used for the burst oscillation. We compare NS spins of both 
300 Hz (Fig. Et) and 600 Hz (Fig. Eb) (dashed lines). 
Generally, a faster spin shows larger frequency drifts for these 
modes because A oc q 2 . This is not the case for the / = 2, m = 1 
r-mode we favor as the burst oscillation, which shows no de- 
pendence on f2 (thick solid lines). When the cooling is fastest 
at early times, the frequency shifts are dramatic, ranging any- 
where from w 15 Hz (for the Kelvin modes) tow 100 Hz (for 
the rotationally modified g-modes). Such large shifts have 
not been seen in previous frequency searches of X-ray burst 
lightcurves, but may be lurking, undetected. On the other 
hand, because of the late time frequency stability afforded by 
the crustal mode, one could argue that if these modes existed 
they would have already been detected when the frequency 
evolution is weak. The resolution of this question depends on 
the actual oscillation amplitude. If the burst oscillations are 
any indication, the amplitudes are largest soon after the burst 
peak when the frequency shift is greatest. The stable part of 
the frequency might then have an amplitude below detectable 
levels. It remains an important exercise to search the exten- 
sive database of burst lightcurves for signatures of oscillations 
similar to what we describe here. 

7. CONCLUSIONS AND DISCUSSION 

We considered nonradial surface oscillations as a possible 
explanation for the oscillations observed during type I X-ray 
bursts on non-pulsar accreting NSs. In studying the time evo- 
lution of a shallow surface wave, we found that this mode 
changes into a crustal interface wave as the surface cools, 
a new and previously unexplored result. This phenomenon 
allows us to match both the observed frequency shifts and 
the high stability of the asymptotic frequency of burst os- 
cillations. To find frequencies in the observed range, the 
modes must be highly modified by rotation, resulting in the 
NS spin occurring merely w 4 Hz above the burst oscill ation 
frequency. Following similar reasoning to iHevI ( 120041) . we 
find that the majority of the rotationally modified modes can- 
not match all of the observed properties of burst oscillations, 
because they predict too large of drifts and/or drifts of the 
wrong direction. The I = 2, m = 1 buoyant r-mode (A w 0.1 1) 
overcomes all of these difficulties, and for this reason it is the 
favored mode to explain the burst oscillations. Perhaps not 
coincidentally, this mode has other attractive qualities such 
as a wide latitudinal eigenfunction that may make it easier to 
observe. 

An important mystery that we did not resolve here is why 
burst oscillations in the non-pulsars are only seen in the short 
~ 2 — 10 s bursts (not counting the superbursts), a decay 
timescale indicative of He-rich fuel at ignition. While Models 
1 and 2 show frequency and amplitude evolution similar to 
what has been observed for burst oscillations, Model 3 does 
not. This is not unexpected because the first two (especially 
Model 1) reflect what we expect for He-rich bursts, but it is 
at the same time perplexing because we have no reason a pri- 
ori to explain why burst oscillations as from Model 3 are not 
observed. There must then be some phy sical mechanism for 
why th ey are not seen. In the study by Cumming & Bildsten 
(2000), they find that modulations in the flux at the burning 
depth can be washed out if the thermal time is too long in 
the hydrostatically expanding and contracting surface layers. 
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FIG. 12. — Frequency evolution of the lowest angular-order, rotationally modified modes from Model 1. All modes are 1 = 1, unless otherwise noted, where 
/ refers to the latitudinal quantum number in the slowly spinning limit. In every case the radial structure is the single-node mode as we focus on for the burst 
oscillation. The mode frequencies are denoted by solid lines, with the NS spin as a dashed line. In (a) we consider a spin of 300 Hz and find that its frequency 
shifts are somewhat smaller than for a 600 Hz spin as shown in (b). Even so, all the frequency shifts are much larger than the / = 2, m = 1 r-mode that we favor 
as the burst oscillations (which is also plotted for comparison; thick solid line). 



This effect is dependent on the mean molecular weight so that 
it is weak for He-rich bursts, but fairly strong for mixed H/He 
bursts. It is possible that a similar phenomena occurs just 
above the nonradial modes that we study here. 

There remain many properties of burst oscillations that are 
not stud ied in our work here (for example, see t he observa- 
tions of M uno. Ozel. & Chakrabartvl l2002l |2003 ). The char- 
acteristics of the oscillating light curves, the energy depen- 
dence of the pulsed amplitudes, and the phase shifts of differ- 
ent energy components all can be used to learn about the ac- 
creting NSs that are exhibiting these burst o scillations. Such 
investigations have just b een initiated by iHevll J2005I) and 
iLee &Strohmaverl (120051) . 

The most exciting implication of our study is that the 
asymptotic frequencies and drifts are directly related to prop- 
erties of the surface NS layers. The crustal interface wave 
depends on the attributes of the crust, which in turn is a prod- 
uct of the X-ray burst properties, and especially its superburst 
properties. Superbursts are deep enough to directly impact 
the composition of the NS ocean, and eventually, as material 
advects down with the accretion flow, the crust. The accret- 
ing low mass X-ray binary 4U 1636-536 has exhibite d oscil - 
lations during a superburst (Strohmaver & Markwardt 2002), 
and during X-ray bursts both before and after (Zhang et al. 
1997; M. Muno 2004, private communication), all with the 
same asymptotic frequency. Energetically, superbursts could 



easily excite crustal waves since they are roughly a thousand 
times more energetic than normal X-ray bursts. For the fre- 
quency to be the same in each occurrence the crustal proper- 
ties could not have changed during this time. This is possible 
as long as 4U 1636 — 53's superburst recurrence time is short 
enough that the crust keeps a steady state composition. A 
crust at « 10 13 g cm" 2 requires «s 7 yrs at an accretion rate of 
10~ 8 Mq yr" 1 to be replaced by new, accreted material. Such 
a timescale favors a crust that has not changed significantly 
given that 4 U 1636-53 has pro duced three superburst s in the 
last 4.7 yrs JWiinandJ200lHStrohmaver & Markwardl2002l 
iKuulkers et al.l2004l) . 
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Note added in proof. — It has been brought to our attention by Dong Lai that our arguments in SI4.3I for adiabatic evolution 
through the avoided crossing may not be correct. We currently use the condition that adiabatic evolution occurs as long as the 
mode frequency is much larger than the cooling rate. If the avoided crossing instead acts analogously to neutrino oscillations (J. 
Bahcall 1989, Neutrino Astrophysics [Cambridge: Cambridge Univ. Press]) or photon propagation on magnetized neutron stars 
(D. Lai & W. C. G. Ho 2002, ApJ, 566, 373) the correct comparison would be the frequency difference at the avoided crossing 
versus the cooling rate. If this is indeed the case, then the surface mode will not solely couple to the crustal interface mode. 
Instead, some energy would remain in the surface mode as the frequencies cross. We are currently studying this problem. 
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Properties of X-ray Burst Oscillations from Non-pulsars 
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"Asymptotic frequency. 

b The range of frequency drifts seen from each object. 

c Estimate of the decay time of cooling in those bursts that exhibit oscillations. This should be considered 
as only a rough estimate since different authors use different methods of measuring this quantity. The main 
conclusion to be inferred here is that these are all short, He-like bursts. 

''Whether or not radius expansions are observed in those bursts that exhibit oscillations. Aql X- 
1 usually shows radius expansions with burst oscillations, but there is one case when it did not (see 
Muno . Galloway. & Chakraba rty|2004l) . 

e The persistent luminosity as a percentage of the Eddington luminosity. Each range of measurements is 
shown for those objects with multiple observations. We use Lejj = 2.5 X 10 38 erg s~' to be consistent with 
IFordetal]<200TI) . 

f This drift could be interpreted as preceding the burst peak. 



TABLE 2 
Cooling Neutron Star Models 
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